arXiv:1506.03450vl [astro-ph.SR] 10Jun2015 


Draft version June 12, 2015 

Preprint typeset using ET^X style emulateapj v. 08/22/09 


IMPRINTS OF EXPANSION ONTO THE LOCAL ANISOTROPY OF SOLAR WIND TURBULENCE 

Andrea Verdini^ 

Dipartimento di fisica e astronomia, Universita di Firenze, Firenze, Italy. 


Roland Grappin 

LPP, Ecole Polytechnique, Palaiseau, France. 

(Dated: June 12, 2015) 

Draft version June 12, 2015 

ABSTRACT 

We study the anisotropy of Il-order structure functions defined in a frame attached to the local mean field in 
three-dimensional (3D) direct numerical simulations of magnetohydrodynamic turbulence, including or not the 
solar wind expansion. We simulate spacecraft flybys through the numerical domain by taking increments along 
the radial (wind) direction that forms an angle of 45^ with the ambient magnetic field. We find that only when 
expansion is taken into account, do the synthetic observations match the 3D anisotropy observed in the solar 
wind, including the change of anisotropy with scales. Our simulations also show that the anisotropy changes 
dramatically when considering increments oblique to the radial directions. Both results can be understood by 
noting that expansion reduces the radial component of the magnetic field at all scales, thus confining fluctua¬ 
tions in the plane perpendicular to the radial. Expansion is thus shown to affect not only the (global) spectral 
anisotropy, but also the local anisotropy of second-order structure functions by influencing the distribution of 
the local mean field, which enters this higher-order statistics. 

Subject headings: The Sun, Solar wind. Magnetohydrodynamics (MHD), Plasma, Turbulence. 


1. INTRODUCTION 

The solar wi nd is known to be a turbulent medium since 
many decades (lColemanlll968l) and is probably the best ex- 
ample of natural turbul ent laboratory in astrophysics (e.g. 
iBruno & Carbond 120131) . Turbulence shows most of the 
time a non zero global mean field which should lead 
to an anisotropi c cascade with the spectrum being axisym- 
metri c around it (IMontgomerv & TurnedI 19811: IShebalin et al.l 
119831: iGrappinl Il986l) . As the angle between and the 
radial direction varies, a spacecraft embedded in the radial 
solar wind samples data in different directions with respect 
to the mean field. This allows one to measure the correla¬ 
tion function in two dim ensions, which has the chacteristic 
of an anisotropic cascade (iMatthaeus et al.ll 19901: iBieber et^ 

rr^lfj^et alJl2nn5tlHamilton et al Jl2nn8h. 

However, the axisymmetry assum ption has been found 
to break down i n several works (ISaur & Biebed 119991: 
iNarita et ^120101: IChen et ani2012l) . This may result from 
(a), considering scales large enough for the expansion to play 
a role and/or (b), considering anisotropy with respect to the lo¬ 
cal mean field instead of the global mean field. While having 
a reference frame attached to the glob al is preferable for 
studying the turbulent dissipation (e.g. lVerdini et al .1120151) . a 
reference frame attached to the local, scale-dependent, mean 
field allows one to reveal the effect of local dynam¬ 
ics in magnetohydrodynamic (MHD) turbulence. In the lat¬ 
ter case a different scaling in the two directions parallel 
and perpendicular to was found both in direct numerical 
simulations (DNS) (e.g. lOio & Vishniacll200()l: iMilano et al.l 
120011: iBeresnyak & Lazarianl120091: iGraopin et al.ll2013l) and 


in the sola r wind (e.g. iHorburv et alJ 12008 : IPodestal 
iLuo & Wul [Mol: IWicks et al.l iMol iMll 1201 21 


2009; 


2013 


^ Solar-Terrestrial Center of Excellence - SIDC, Royal Observatory of Bel¬ 
gium, Bruxelles, Belgium. 


IChen et all IIMiltlHe et aillMnh . 

Deviations from axisymmetry (in the form of three distinct 
scaling laws) appear when consideri ng two perpend icular di¬ 
rections instead of a s ingle one (see lBoldvrevll206^ . In their 
recent measurements ICtien et (1201 2b show how the small 
scale ordering of the structure functions (SF) is completely 
modified in the solar wind when passing from small to large 
scales. While the small-scale anisotropy is roughly compati¬ 
ble with t hree-dimension al anisotropic phenomenology of tur¬ 
bulence (lBoldvrevll2006b . the large-scale anisotropy has no 
explanation so far. 

In this Letter we focus on the large-scale ordering and 
we explain it with phenomenological arguments supported 
by DNS of MHD equations modifie d to include expan¬ 
sion (expanding box model or EBM, iGrappin et aki 1199^: 
iGrappin & Vellil T ~9^ . The EBM has been recently used 
(iDoing et al .1120 14b to show the scale-dependent competition 
between two axes of symmetry, the mean field axis and the 
radial axis. Here we show that it is able to reproduce both 
the large and small scale anisotropy of the SF along the three 
ortogonal direction defining the frame attached to the local, 
scale-dependent mean field. 

2. SIMULATIONS AND PARAMETERS 

We follow the evolution of a plasma volume embedded in 
a mean flow with constant speed. Turbulent evolution with 
distance is thus modeled as decaying, unforced turbulence. 
Two runs are analyzed: run A assumes a uniform parallel 
mean flow, run B assumes a radial mean flow, as is the so¬ 
lar wind. The full MHD equations (continuity, induction, ve¬ 
locity and energy equations), are integrated in time with a 
pseudo-spectral code on a grid of 1024^ points. For run B, 
the MHD equations are modi fied to incorporate expansion, 
becoming the EBM equations (iDong et al .1120 14b . In the fol¬ 
lowing, velocities are normalized to the initial rms amplitude 
of velocity fiuctuations Urms^ lengths to the box size L^, and 




























































2 


(a) Global frame (b) Local Frame 

(x,y,z)^(R,T,N) 



Fig. 1. — Panel (a). Reference frame for simulations at f = f. The R, T, N 
directions correspond to the x, y, z coordinates of the numerical domain. 
The mean magnetic field Bq forms an angle Obo = 45^ with the direction 
of increments i = £x connecting the fiuctuations Bi and B 2 . Panel (b). 
Local reference frame for computing the 3D anisotropy of structure functions. 
The axes A and the angles 6b, 6sb± are defined for each couple of 
points by the local magnetic field Be = (Bi + B 2 )I 2 and by the fiuctuation 
6B = B 2 — Bi. 


time to the initial eddy turnover time jlnurms • The 

magnetic field B is also expressed in unit of Alfven speed, 
B! A/47rpo, with po being the average density. 

We first define the expanding run B. The reference frame, 
X, y, z, is aligned with the R, T, N coordinates of the he¬ 
liocentric reference frame (Figure [T^). The domain is ad- 
vected by the solar wind at a constant speed Vsw' dif¬ 
ferent times t correspond to different heliocentric distances, 
R(t) = Vswt, with RP being the initial position of the 
simulation domain. During advection, the domain inflates 
anisotropically: the radial dimension does not change with 
time, while the lateral dimensions scale Ly^L^ oc R{t). 
The rate of inflation is set by the expansion parameter, e = 
= (LV2nur^,)l(Ro/Vsw), where = R^Ysw is the 
initial expansion time. We fix the initial heliocentric distance 
R^ = 0.2 AU, the lateral dimension of the numerical domain 
Ly = = LP = R^/5, and the ratio UrmsIVsw = 1/4;:, yield¬ 

ing finally e = 0.4. The initial aspect ratio is R^ = Ly/L^ = 
1 / 5 so that at R = 1 AU we have a cubic numerical domain 
= R^ = 5LP = 0.2 AU. The conservation of magnetic 
flux implies Bq oc (1/R, 1,1), we thus impose an oblique ini¬ 
tial mean field, Bq = [1,1/5,0], to have an average Parker 
spiral angle of 45^ at 1 AU. Finally we set equal viscosity, 
resistivity, and conductivity, y = rj = k = 6 10 “^ and allow 
the coeffcients to vary as l/R to cope with the damping of 
fiuctuations due to expansion. 

For the non-expanding run A (6 = 0), we choose = Ly = 

L^^=R^,Bo = [V2, V2,0],y = = 1.1 IQ-P 

The fiuctuations B and u are initialized in the same way in 
both runs, as a superposition of modes with random phases, 
with the velocity being divergence-less. Their spectra fol¬ 
low a bi-normal distribution in the Fourier space, of widths 
cr_y = 4k^ and cr\\ = (r^jA for wavevectors perpendicular and 
parallel to Bq respectively (k^ = InILP). The initial eddy- 
turnover time is thus four times smaller in the perpendicu¬ 
lar directions than in the parallel direction. In the expanding 
case, this reduces the expansion effects in the directions per- 
pedicular to the radial. The magnetic and kinetic fields are at 
equipartition, Brms = i^rms = 1, subsonic (the sound speed is 
Cs ~ 7), and have statistically vanishing correlation (m -.B) ~ 0 
(no imbalance between the Elsasser modes). 


To compute the anis otropy of structure functions we use the 
procedure described in IChen et al.l (120121) . For each couple of 
points JCi, X 2 separated by the increment i = JC 2 - JCi, we 
define the local mean field as Be = ll2(Bi B 2 ) and the fluc¬ 
tuating field as SB = B 2 - Bi, where ^ 1,2 = ^(-^ 1 , 2 ). The lo¬ 
cal scale-dependent reference frame, shown in Figure[TJ), has 
the vertical axis f oriented along Be, the first perpendicular 
axis ^ oriented along the perpendicular fiuctuation direction 
SBj_ oc Be X [SB x Be], and the second perpendicular axis A 
perpendicular to both Be and SBj_ . In this reference frame, the 
polar and azimuthal angles 6b and Osbi. define the direction of 
increment with respect to the local mean field. For each pair 
of points the 5-trace structure function, SF = \SB\^, is accu¬ 
mulated in 5^ bins for 6b, 6sb± ^ [0^, 90^], and then averaged 
in each bin (we reflected below 90^ any angles larger than 
90^). Increments, except when otherwise stated, are com¬ 
puted along the x direction i - corresponding to space¬ 
craft flybys along the radial direction in the solar wind frame, 
as in observations. 

We first present the results of the flyby analysis on simu¬ 
lated data at f = 2.8 for run A and at f = 10 for run B. We 
then show how the anisotropy evolves in time in the two runs. 
While f in run B is chosen to reproduce data at R = 1 AU, the 
choice for run A is arbitraty. Homogeneous runs evolve more 
rapidly than expanding runs since, in the latter, fiuctuations 
are damped by both turbulence and expansion and so the non¬ 
linear time increases more rapidly. We thus chose a different 
time in run A, after the peak of current density {t ~ 2.4) but 
not too late, in order to have a Reynolds number Re 1200 
and BrmsIBo ~ 0.9 comparable to those of run B at f = 10 
(Brms I Bo -1.5 andR^ 1300) 

3. RESULTS 

In Figure [ 2 ^ we plot the SF of the homogenous run A as a 
function of wavenumber k = 2nIix for three couples of angle 
6b, 6sbx corresponding to the directions A, f in the local 
reference frame of Figure [TJd. Increments are taken along the 
X direction, which forms an angle of 45^ with the mean field 
Bq. At large scales k < 8 the SFs have comparable energy in 
the three directions. At small scales, 8 < k < 60, we have 
SF(/1) > SF(^) > SF(f) with the the following approximate 
scaling SF oc A^^^, (the power-law range is actually 

smaller in ^ and f). 

In panel (b) we show the same plot for the expanding run 
B. Its overall structure differs completely from run A. Now 
SF(/1) and SF(^) have parallel profiles roughly proportional to 
/fi/^, in the small-scale range, 10 < k < 50. SF(/I) is 
dominant everywhere, while SF(f) passes from almost domi¬ 
nant at large scales (k < 10) to subdominant at small scales, 
where the ordering is the same as for the homogeneous run A. 

Following IChen et aT] (120121) . another viewpoint of the 
anisotropy in the expanding case is given in Figure O where 
we plot the isosurfaces of constant SF power at three different 
levels (marked as dashed lines in Figure [Jb), corresponding, 
from left to right, to smaller and smaller scales. For a given 
value of the isosurface, its shape indicates the correlation of 
fiuctuations along the three directions of the local frame, and 
can be roughly thought as a statistical eddy shape. Since 
SF = \SB\^ measures the power in the anticorrelation, the 

^ The Reynolds number is computed as Re = {LinjlLsissP^^, where Linj = 
(Sn/AE) 2k ^(k)/k and Ldiss = D = 2k y^^E(k) is the dissipation 

per unit mass and E(k) is the omnidirectional spectrum. 














3 


(a) Run A, £ = 0 


(b) Run B, e = 0.4 
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Fig. 2.— Second order SF for the non expanding run, panel (a), and expanding run, panel (b) (compare to Figure 1 in IChen et al.ll2012l) . as a function of 
k = 2nllx- SF are accumulated at three different couples of angles, 6b and 0 , 55 ^, corresponding to the parallel direction (^, red line) and the two perpendicular 
directions (A in blue, ^ in green) respectively (see Figure [TJ)). The dotted lines are a reference for the scaling k~^^^ and k~^ in panel (a) and for k~^ 

in panel (b). The horizontal dashed lines in panel (b) mark the levels at which isosurfaces are drawn in Figure [3] 






Fig. 3.— Run B. Isosurface of constant SF power at three levels, logSF = -0.55, - 0.85, - 1.4, corresponding to smaller a nd smaller scales r espectively 
(horizontal dashed lines in Figure|2). Colors indicate distance from the origin to help the 3D visualization. Compare to Figure 4 in IChen etdl <20121) . 





more energetic is the SF along a given direction, the smaller 
its correlation. In Figure [3l the smallest correlation (smallest 
elongation of the isosurface) is always in the A direction, but 
the direction of the largest correlation changes with scales. At 
large scales (Figure [3^), the eddy is more elongated in the ^ 
direction, corresponding to SBj_, at small scales (Figure[3]:) it 
becomes more elongated along the f direction, corresponding 
to Be. 

The anisotropy shown in Figure IJJ) and Figure [3] are in ver y 
good agreement with the observations of lChen et al.l (l2012h . 
In Figure |2t) the small-scale anisotropy (10 < k < 50) is 
roughly consistent with critical balance, while the large-scale 
behavior (k < 10) requires some more explanation. 

The large-scales ordering of SFs is related to the compo¬ 
nent anisotropy of the magnetic fluctuations that originates 
from the selective damping induced by the expansion. The 
component anisotropy is shown in Figure [4^, where we plot 
the reduced energy spectra of the x, y, and z com¬ 

ponents of the magnetic fluctuations compensated by for 
runs A and B. While in the non-expanding case energy is dis¬ 
tributed isotropically among the components, in the expand¬ 
ing case the radial component is at least a factor 2 smaller. 


at all scales. This behav ior is consistent with observations 
(iHorburv & Baloghll2001[) and is generally found in expand- 
ing runs (see iDong et al.ll2014h . The link between the compo¬ 
nent anisotropy and the SF anisotropy is shown in Figure [4j) 
where we plot their evolution with time. The SF anisotropy 
is quantified as the ratio SF(0/SF(f) at k = Inj lx = 8 (see 
lEsquivel & Lazarianll201 lb iBurkhart et al.ll2014l for a similar 
analysis on global SF), while the component anisotropy is 

evaluated as A = ^(E"^ -f Ey)l2E^ at k^^ = 8. For the ho¬ 
mogenous run A, both ratios are about constant and close to 
the value of 1 (isotropy), while in run B both ratios increase 
steadily and approximately with the same rate. Thus, in the 
expanding case as the heliocentric distance increases the mag¬ 
netic fluctuations are more and more conflned in the y, z plane 
(the T,N plane), and so also the local mean held will prefer¬ 
entially lie in this plane. 

We now show that when the SF is sampled along the ra¬ 
dial direction, the above component anisotropy, E^^,E^^>E^, 
implies that the SF has a different power along the three direc¬ 
tions A, f defining the local reference frame. Consider two 
vectors Bi, Bj at positions xi and X 2 and indicate with B^^ 
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(a) Energy spectra (b) Anisotropy 




Fig. 4.— (a) Energy spectra at f along the radial direction k;^ for run A (top curves) and run B (bottom curves) for the three R, T, N components of the magnetic 
field, Bx, By, Bz (dashed, solid, and dotted lines respectively^ (b) Evolution with time of the structure function anisotropy SE(^)/SF(^)| (thick line) computed at 

k = 2nl£x = ^ and of the component anisotropy A = + E^)/2E^ (thin line) computed at k;^ = 8 for run A (top) and run B (bot). Asterisks indicate the time 

f at which the flyby analysis is presented. 


their components in the y, z plane with a the 

angle between them, and with ^ their projection along x. 
Assume also for simplicity 

B\2~0(2)«Bf =0(l)=l. (1) 

The local mean field and the fiuctuation in the y, z plane are 
given by 


Bf= , 

J\ Bf ^Bf cos Of, 

(2) 

SB-^ = 

A Bf^ — 2Bf cos a 

(3) 


should not be confused with SBj_ that defines the local 
reference frame). These equations simply state that when Bf 
and ^2 aligned (anti-aligned) the local field Bj is large 
(small) and the fluctuating field SB^ is small (large). The ori¬ 
entation of the local reference frame with respect to the fixed 
radial direction x determines which local SF we are measur¬ 
ing: we cumulate the power in SF(/f) or in SF(f) or in SF(^) 
when T or f or ^ lies along x. We thus estimate the leading 
order contributions associated to each of them by considering 
the power \SB\^ associated to the above three orientations: 

• To have power in SF(/f) one needs both B^ and SB to 
lie in the y, z plane, that is, Bf » b'I and 6B^ » 6BK 
This condition is readily satisfied from Equation O, 
provided that Bf « Bj, or Bf ~ Bf and 60"^ < 
a < 120^ in Equations o-©, yielding in both cases 
SF(d) ~ |(5B-^|2 -0(1). 

• The contribution to SF(f) is obtained when B^ is 
aligned along x, that Bf « Bf This happens only 
when Bf ~ Bf and 120^ « a < 180"^ in Equa¬ 
tion (|2]), implying from Equation dS]) a contribution 
SF(0 - \SB^? -0(1). 

• Finally, to have power in SF(^), Bf must belong the 
y, z plane and SB must have a non negligible x com¬ 
ponent SB^^ > SB^I (our minimum angular resolu¬ 
tion is 5^). The conditions SB^ ~ SB^^ ~ (9(2) and 


b\ « Bf ~ (9(1) are satisfied in Equations (O-© 
only when Bf ~ Bf and a « 60^, which yield a power 
SF(^) ~ \SB^\^ ~ ~ (9(2). 

Thus, the geometrical constraint imposed by the component 
anisotropy induced by expansion, favors a local 

anisotropy with SF(/f) > SF(f) > SF(^), as is indeed found at 
large scales in Figure [Jb and in the solar wind observations. 

We finally show how the anisotropy of expanding turbu¬ 
lence changes if one samples increments in directions per¬ 
pendicular to the radial. In Figure [5] we plot the SFs of run 
B computed along Sx and iy, corresponding to the R and T di¬ 
rections, respectively. The SFs are now compensated by 
to highlight inertial-range scales and the corresponding spec¬ 
tral index. Independently of the direction of increments, the 
inertial range extends to smaller scales in the second perpen¬ 
dicular direction (T) than in the first perpendicular direction 
(^). The parallel direction (f) does not show any convinc¬ 
ing scaling, although having a steeper spectrum. The direc¬ 
tion of increments affect the overall anisotropy. In fact, the 
large-scale ordering, characteristic of expansion in panel (a), 
disappears when increments are along the transverse direction 
(panel (b)), the SF becoming basically isotropic for ^ < 10. 
This can be interpreted as a reduced effect of the component 
anisotropy when increments are along the transverse direction 
(a similar behavior is seen for S = not shown). Moreover 
the two perpendicular SFs (blue and green curve) exhibit a 
different spectral index, passing from a slope 1 /2 in panel 
(a) to steeper spectral index 2/3 in panel (b). Note finally 
that the eddy shape for transverse increments is at all scales 
qualitatively similar to the anisotropy in Figure [3]:, and thus 
differs completely from the cases of radial increments shown 
in Figure [3^,b. 

4. DISCUSSION 

We computed the 3D anisotropy of structure functions with 
respect to the local mean field in DNS of MHD turbulence, 
including or not expansion. For homogenous turbulence, the 
SF is roughly isotropic at large scales and develops scale- 
dependent anisotropy at small scales due to the different scal¬ 
ing along the different local directions. The corresponding 
spectral indices are roughly consistent with -1/2, -2/3, and 
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Fig. 5.— SF of expanding run B compensated by as a function of increments computed along the radial directions k = 2nlix (left), and along the transverse 
direction k = 2nliy (right). The two dashed lines are a reference for the scaling k~^^^ and . 


(a) run B, increments along 


(b) run B, increments along 
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-1 in the A, ^ and ^ direc tions resp e ctively . Such an order¬ 
ing has been predicted by iBoldvrevI (120061) and implies that 
SF(/1) > SF(^) > SF(f ) at small scales. The homogenous run 
qualitatively agrees with the above power anisotropy, which 
is stable although the precise slopes vary with time and with 
sampling direction. 

When expansion is taken into account the anisotropy is de¬ 
termined by a well-defined large-scale power anisotropy and 
by a different scaling along the parallel direction (f) and the 
two perpendicular directions (A, ^), the latter now having ap¬ 
proximately the same spectral index -1/2). We found 
that the overall SF anisotropy is a consequence of the com¬ 
ponent anisotropy induced by expansion and that it shows up 
only when increments are computed along the radial direc¬ 
tion. When increments are along the transverse direction, 
the large-scale anisotropy disappears, the eddy shape does 
not change with scales (although the anisotropy increases at 
smaller and smaller scales), and the SFs exhibit steeper spec¬ 
tral scaling. Thus, the measured anisotropy of solar wind tur¬ 
bulence would change for increments in directions other than 
the radial, a situation that may become possible to test with 
Solar Probe Plus in its near-sun orbital phase. 

Let us compar e our resul t s with the observations in the 
fast solar wind (IChen et al.l l2012h . In run B the choice 
of the solar wind speed yields, via the Taylor hypothe¬ 
sis, the spacecraft frequencies corresponding to the radial 
increments in Figure [2] and Figure O For a fast wind 
Vsw ~ 800 km s~\ one gets = Vsw^x/^^ ^ 
and fmax = 512 X f^m -1.5 10“^ Hz. Our initial fluctuations 
have vanishing u, B correlations for easier comparison 
between the expanding and homogenous runs, while fast 
wind has high correlations. Contrary to the non-expanding 
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